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ABSTRACT 

We study a class of nonlocal, but causal, covariant and conserved field equa- 
tions for the metric. Although nonlocal, these equations do not seem to 
possess extra graviton solutions in weak field perturbation theory. Indeed, 
the equations reduce to those of general relativity when the Ricci scalar van- 
ishes throughout spacetime. When a static matter source is present we show 
how these equations can be adjusted to reproduce Milgrom's Modified New- 
tonian Dynamics in the weak field regime, while reducing to general relativity 
for strong fields. We compute the angular deflection of light in the weak field 
regime and demonstrate that it is the same as for general relativity, resulting 
in far too little lensing if no dark matter is present. We also study the field 
equations for a general Robert son- Walker geometry. An interesting feature 
of our equations is that they become conformally invariant in the MOND 
limit. 
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1 Introduction 



Milgrom's Modified Newtonian Dynamics (MOND) is an empirical alter- 
native to invoking dark matter as an explanation for the motions of cosmolog- 
ical systems which are subject to very small accelerations. In one formulation 

— # 

the force FNewt on a point mass m is given by Newtonian gravity, but the 
second law becomes nonlinear in the acceleration a, 

-FNewt = m/if— )a where /i(x) = < 1 ^ X ^ ] ' (1) 
\a / he Vi«l. 

The numerical value of a® has been determined by fitting to the rotation 
curves of nine well-measured galaxies |2J, 

a = (1.20 ± .27) x lCT 10 m/s 2 . (2) 

Because the observed stars and gas contain constituents subject to much 
larger internal accelerations 1 it is preferable to view MOND as a modifica- 
tion of the gravitational force at low accelerations. That is, the second law 
takes its traditional form -Fmond = ma, but the actual force is the following 
nonlinear function of the force predicted by Newtonian gravity, 

^mond = /f^^'l^Newt where f(x) = \ 1 L X > ] ' (3) 

One consequence of MOND is that the velocities of particles in progres- 
sively more distant circular orbits approach a constant Voo which depends 
only upon a galaxy's mass |3 3 . For large radius r the Newtonian force due 
to a galaxy of mass M is FN ew t = GMm/r 2 . Since F^ ewt /m must even- 
tually fall below ao, relation (JHJ) implies that the asymptotic force law is 
-Fmond = my/a GM jr. Applying the second law gives, 



V <*f M = V -f =► a GM = vt o . (4) 

With no dark matter a galaxy's luminosity L should be a constant times its 
mass, although the constant will depend upon the type of galaxy. Thus one 
expects L ~ v^, which is the observed Tully-Fisher relation jl]. 

1 For example, were we to think of neutral Hydrogen as a classical planetary system, 
the proton would be subject to an acceleration of about 10 19 m/s 2 about the atomic 
barycenter! 
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MOND provides an excellent fit to the rotation curves of all known types 
of galaxies, using only the measured distributions of gas and stars and fitted 
mass-to-luminosity ratios for gas and stars. The recent review by Sanders 
and McGaugh [S] summarizes the increasingly compelling data and lists the 
primary sources. It is significant that the fitted mass-to-luminosity ratios are 
not unreasonable [U]. Especially significant is that MOND agrees in detail 
even with low surface brightness galaxies [7||B], objects for which the MOND 
regime (a < ao) applies throughout and for which no detailed measurements 
had been made when MOND was proposed. On the other hand, some dark 
matter must be invoked to explain the temperature and density profiles of 
large galaxy clusters 0. (It it has been suggested that eV^-mass neutrinos 
might provide this without affecting the galactic results |1U|.) And a very re- 
cent analysis of data from the Sloan Digital Sky Survey asserts that satellites 
of isolated galaxies violate MOND when care is taken to exclude interlopers 

The successes of MOND might be a numerical coincidence ^2] • They may 
result instead from galaxy formation and evolution flowing towards some yet- 
unrecognized attractor solution through conventional physics. Or they might 
signal a real modification of gravity in the regime of very low acceleration. 
Although the issue cannot be decided at this time, the number and quality of 
relevant observations which are technically feasible indicate that it may 
be resolved in the near future. If MOND could be convincingly embedded in 
a larger, metric formulation of gravity, that theory could already be tested 
against a wide variety of data from lensing [T3J [T5|, cosmology [T£|, and 
structure formation ^7j. As it is, one must make assumptions, whereupon 
comparison with observation is as much a test of these assumptions as it is 
of MOND. It is therefore opportune to consider what sort of fundamental 
theory might reduce to MOND in the appropriate limit. 

There is a satisfactory nonrelativistic potential formulation for the rela- 
tion between the MOND force Fmond = — V0 and the mass density p m . It 
is given by the Lagrangian of Bekenstein and Milgrom ^H] , 

£ = -Prn<p-^F(^^) where p(x) = F'(x 2 ) . (5) 

This is very important because it establishes that MOND conserves energy, 
momentum and angular momentum. However, it does not extend the theory 
sufficiently to test lensing and cosmology. A generally coordinate invariant, 
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scalar-metric extension exists but it contains dynamical scalar degrees of 
freedom which can propagate acausally [T8] . 

We wish here to consider a different possibility. Suppose that general 
relativity really is the fundamental theory of gravity, but that its effective 
action contains large quantum corrections from infrared virtual particles. 
Weinberg showed that infrared effects in quantum gravity are no stronger 
than those of QED for the case of zero cosmological constant ^H]- However, 
a nonzero cosmological constant would preserve the graviton's masslessness 
while subjecting it to interactions of canonical dimension three. Recall that 
infrared effects become stronger as massless particles are coupled with lower 
dimension interactions, and that they are already nonperturbatively strong 
for massless gluons with the dimension four coupling of QCD. 

These considerations are the basis for a daring proposal to simultaneously 
resolve the (old) cosmological constant problem (213 E] and provide a natural 
model of inflation in which scalars play no part. The idea [22J is that the bare 
cosmological constant is actually GUT-scale, which leads to an initial period 
of inflation during the early universe. What brings inflation to an end is 
the gravitational attraction between the ever-increasing numbers of infrared 
gravitons ripped out of the vacuum by the rapid expansion. One can follow 
this process as long as the slowing is weak, and explicit computations confirm 
that the effect must eventually become nonperturbatively strong [23J. 

If this proposal is correct, the post-inflationary universe would be de- 
scribed by the nonlocal effective action which prevails after the breakdown 
of perturbation theory. We cannot compute reliably in this regime, but noth- 
ing prevents one from making guesses about the form of this effective action. 
The simplest class of guesses which give a plausible end for inflation involve 
acting the inverse covariant d'Alembertian on the Ricci scalar [23]. We will 
refer to this as the small potential, 



(We use a spacelike metric with R, u = T^ p - T»^ v + T>> p(T T% - T^T^.) 



This paper will not consider inflation or the cosmological constant. We shall 
rather explore embedding MOND in a nonlocal Lagrangian of the form, 



A desirable feature of this class of models is that it involves only the 
metric. Although the field equations are nonlocal, they do not seem to possess 



(p[g] = — R where 
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additional graviton solutions in weak field perturbation theory. To see this, 
expand the metric about flat space as usual, 

9tw = Vtiv + hpv . (8) 

It is an elementary exercise to show that the Ricci scalar is, 

R=h^-hYu + 0(h 2 ), (9) 

where graviton indices are raised and lowered using the Lorentz metric as 
usual. Now impose de Donder gauge, 

h\» ~ \h^ v = , (10) 
to show that the small potential is local in the weak field limit, 

<p[v + h] = -±h", l + 0(h 2 ). (11) 

Since the Lagrangian depends upon the first derivative of (p, this theory 
cannot acquire higher derivative solutions in weak field perturbation theory. 
Indeed, all solutions to the source-free Einstein equations must be solutions 
to this theory because they have R = throughout spacetime, which implies 
if = as well. The correction term in (J7J) changes only gravity's response 
to sources, without adding new weak field dynamical degrees of freedom. So 
these models seem worthy of closer study. 

In section 2 we derive the form a static and spherically symmetric met- 
ric must take to reproduce galactic rotation curves in the MOND regime. 
We also work out the deflection of light in the weak field limit. In section 
3 we employ a formal shortcut to obtain the causal field equations which 
would come from a Schwinger-Keldysh effective action [23]. We also explic- 
itly demonstrate conservation. The nonlocal field equations are specialized 
to a static, spherically symmetric geometry in section 4. We first prove, for 
general interpolating function !F(x), that any MOND terms must drop out of 
the formula for the deflection of light. This is the decisive point of our anal- 
ysis because the existing lensing data can only be made consistent with the 
absence of dark matter if MOND terms contribute ^Hj- No member of this 
class of models can therefore represent a phenomenologically viable extension 
of MOND. However, it is possible to choose the interpolating function T{x) 
to reproduce MOND, and we derive the required form. Section 5 explores 
the model's impact on cosmology, more to illustrate what can be done with 
a complete theory than because this model is viable. In section 6 we discuss 
more general models which might avoid the lensing disaster. 
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2 Phenomenological constraints 



MOND was developed as an alternative to halos of dark matter surrounding 
nearly static, cosmological sources. Far from such a system one would expect 
the asymptotic geometry to be spherically symmetric and static: 

ds 2 = g ilv {x)dx il dx v = -B(r)dt 2 + A(r)dr 2 + rW . (12) 

We first show that rotation curves fix B(r) in the MOND limit but tell us 
nothing about A{r). We then discuss the magnitude of the effect. The 
section closes by working out the angular deflection of light for a general 
class of weak field limits. 

Let x^it) be the worldline of a test particle moving in the geometry 
described by f!12|) . If only the gravitational force is significant the particle's 
worldline obeys the geodesic equation: 

r(t) + T^[ X (t))x p (t)r(t) = o. (is) 

The nonzero connection components derived from (|12|) are, 

r* — ^ r r — ^' r r — ^ r r — r r r — r a ^ 2 (a\ 
r ~ 2B ' ~~ 2A ' r ~ 2A ' 6 ~ ~a > <t><t> ~ T ^ > ' 

rV = \,?*+r = \, = - sin(0) cos(^) , T\ 6 = cot(0) . (14) 
Now specialize to the case of circular motion, 

(V,X r ,X e ,/) = (ct,r,|,#t)) • (15) 

With our connection (|14|) the ji = t and // = 9 components of the geodesic 
equation are tautologies. The \i = <f> component just says <p is constant. Only 
the /i = r component of the geodesic equation is nontrivial, 



B' r 6 2 

2A-17? = - (16) 



Note that A(r) factors out! For circular orbits the velocity has the Eu- 
clidean relation to the angular velocity, v = r<j). In the MOND limit v 2 
approaches the constant v 2 ^ = ^ciqGM, so the MOND limit for B(r) must 
obey, 

B{T) ^lJ^*. (17) 

r V c 4 
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It is well at this point to consider the size of things. A large galaxy might 
have a mass in stars and gas of M ~ 10 11 x M ~ 10 41 kg. Such a galaxy 
would enter the MOND regime at a radius of about, 



tf gal ~ J— ~ 10 20 m. (18) 
V «o 

For weak fields we can write, 

A(r) = l + a(r) , B(r) = 1 + b(r) , (19) 

where \a(r)\ <C 1 and \b(r)\ < 1. A useful phenomenological ansatz for the 
asymptotic behavior of the weak fields is, 



,s x GM a GM GM / a GM 

— ^ + £ iv — > b ^ — + e2 V — ln Vi?g 

(20) 

We have just seen that MOND predicts e<i = 2 but says nothing about the 
other parameters. With just the isolated galaxy, general relativity gives 
5\ = —62 = 2 and e± = €2 = 0. If an isothermal halo of dark matter is added, 
whose density is chosen to reproduce v ^ = \/a GM, general relativity gives 
5\ = —5 2 = ei = £ 2 = 2. 

One might worry that the logarithmic growth of b(r) in ()20j) must even- 
tually invalidate the weak field approximation but this is not a practical 
concern. For the large galaxy considered previously the small parameters 
multiplying the <5's and e's are, 



gm ^ xl0 _ 8 JZ&L-IV*. (21) 

err r V c 4 



The difference b(r) from the onset of MOND all the way to the current 
horizon (i?hor ~ 10 26 m) is, 

b(R ho r) ~ KR&) ~ -5 2 x 10- 6 + e 2 x 10~ 5 . (22) 

The weak field regime is therefore applicable throughout the Hubble volume. 

Another numerical fact worthing noting is that the natural length as- 
sociated with the MOND acceleration is larger than the Hubble radius, 

c 2 

10 27 m . (23) 

a 
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This has the important consequence that ra Q /c 2 is negligible on galaxy and 
galactic cluster scales, so that powers of r do not necessarily distinguish 
"large" and "small" terms. Consider, for example, two possible contributions 
to the weak fields, 



a GM 



In 



GM 

c 2 r 



(24) 



C" 1 VRgal, 

The left hand term is negligible with respect to the right hand side, even 
though the latter falls off with r whereas the former actually grows. 

We can now consider what MOND says about lensing. The angular de- 
flection of light can be expressed in terms of the turning point Ro, 

A?(r) 



Acf) 



dr 



Ro 



7T 



(25) 



gCgo) 
B(r) 



This can be expanded in powers of the weak fields and then simplified using 
the change of variables r = R sec(6 l ), 



A0 



dr 



1 



Ro 



1 



a{r) 



Ro ) 



l &Cgo 
2 



b(r) 



+ 



7T, 



d0|aLRbsec(0)J - csc 2 (#) [b[Ro) -b(Rosec(9) 

Substituting the general ansatz (|2Dj) gives, 

GM 



5l - 6' 



c 2 R ( 



+ 61 



7T 

r3, 2 



a GM 



+ 



(26) 



(27) 



(28) 



Without dark matter, general relativity (5\ 



-6* 



2 and e x = e 2 = 0) 



gives too little deflection at large -Ro to be consistent with the frequency 
of lensing by galaxies. General relativity with an isothermal halo of dark 
matter (6% = —5 2 = t\ = e 2 = 2) is consistent with the existing data [T3] . 
For MOND to be similarly consistent with the data requires the sum (€1 + 62) 
to be positive and of order one. 



3 The nonlocal field equations 

The ostensible purpose of this section is to derive the field equations associ- 
ated with (JZJ). We are not going to quite do that for the very good reason 
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that one does not get causal field equations by varying a temporally nonlocal 
action. To see this, consider the quadratic part of the effective action for a 
real scalar field (f>(x) in flat spacetime, 

' f d 4 y<f>(y)[d 2 - rn 2 ]<j>(y) - l [ d% [ rf 4 #(y)M 2 (x; y)<j>{z) . (29) 



2 J " r W/L JT w/ 2 . 

Here M 2 (y; z) is the scalar self-mass-squared, which is symmetric under in- 
terchange of y^ and z^. Taking the variational derivative and ignoring surface 
terms gives, 

JIM = [d 2 - m 2 ]0(x) - / d 4 x'M 2 (x; x')(j>{x') . (30) 

0(p\X) J 

If M 2 (x;x') for some x ,fl = x M — A M in the past of x^ then it is also 
nonzero for = x^ + A^ 1 in the future of x^. Another embarrassing property 
of the conventional effective action field equations is that their solutions are 
not typically real even if the associated field operators are Hermitian. This 
is because the effective action field equations are obeyed by in-out matrix 
elements of the field operators, and the "in" state may not evolve into the 
same "out" state. 

Both causality and reality can be enforced by employing the field equa- 
tions of the Schwinger-Keldysh effective action [23]. These are obeyed by the 
expectation values of the field operators in the presence of the same state, so 
Hermiticity of the operator implies reality of the solution. In the Schwinger- 
Keldysh formalism there are two background fields: </>+(x) for the evolution 
forward from the initial state and <p-(x) for the reverse evolution back to it. 
The quadratic part of the scalar action has the form, 

4>-\ = \f d 4 yMv)[d 2 - m 2 }^) d'y<P^y)[d 2 - m 2 ]0_(y) 



+ 2J dy J d*z{My)Mt4y;z)Mz) + MyWUy-,z)<i>-(z) 

+ cf>4y)M 2 Jy;z)Mz) + MyW 2 Mz)<t>-{z)} > (31) 



where M+ + (y; z) and M 2 _(y; z) are symmetric, and M+_(y; z) = M 2 + (z; y). 
The various ± self-mass-squared functions are related by simple rules which 
alter the ie terms in propagators and change the sign of some vertices 
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The Schwinger-Keldysh field equations are obtained by varying with respect 
to either <p + or 0_, and then equating the two fields after the variation [25] . 



+1 



[d 2 — m 2 ](f)(x) — 



d 4 x'[Ml + (x;x')+Ml_(x;x')}^(x') . (32) 



0<f>+[x) 

It turns out that M^__(x;x') is exactly equal and opposite to M ++ (x;x') 
whenever x ,fl is not in or on the past light-cone of x^. When x ,fl is in the 
past of x M the two terms are complex conjugates, which ensures reality. 2 

Since we do not actually derive either the conventional effective action 
or the Schwinger-Keldysh version, we will simply employ a trick to extract 
causal field equations from (J7J). For the purposes of this paper one may 
as well regard these equations — rather than ((7|) — as defining the model. 
The trick is to act the nonlocal operators backwards whenever they would 
ordinarily act upon the variation. For example, when f[g] is any functional 
of the metric we write, 



Sy[g] 

$gl»> 



fig] 



i sn 



n Sg» v n 
sn sr 



SR 



nSgt 

-Jig] ■ 



It is useful to recall the standard result for varying the Ricci scalar. 
SR(y) 



Sg^ v (x) 



D^D U 



x 



(33) 
(34) 

(35) 



(Here is the covariant derivative operator.) We recall also the definition 
of the stress-energy tensor from the variation the matter action S m , 3 



T, 



2 5S„ 



-g Sg 



(36) 



Taking IQttGc 4 / \/ — g times the variation of our nonlocal action (JZJ) 
the sense of our trick — gives the following field equations, 

SnGc-%, = 2[$.^ - g^US] + G^l - 2$] 

,2 



111 



+ 



(37) 



2 Full details can be found in a recent paper on the vacuum polarization of scalar QED 
in a locally de Sitter background 26 . 

3 Tliis should be local because only gravitons combine the properties (natural massless- 
ness and the absence of conformal invariance) needed to produce a strong infrared effect 
during inflation |22j . 
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The symbol $ [g] in this and subsequent formulae denotes the large potential, 



(38) 



Because we have not really derived (|37jl it is worth while to explicitly 
verify the important property of conservation. The covariant divergence of 
each term is, 



((^[1-2$] 

(-iM ! 



+2R»% V . (39) 

-2R»% + m tll , (40) 

-<^n$ - m,^ , (4i) 

+y.^y> v T' + <^0$ , (42) 

-VwV'T' , (43) 



and they do sum to zero. The equations (|37j) are manifestly covariant. If 
!/□ denotes the retarded Green's function they are also causal. 



4 Spherically symmetric, static sources 

This is the key section of the paper. We begin by working out the small 
potential (JUJ) and the large potential (|S5|) for a general spherically symmetric 
and static metric ([12)1 . Then we give the two independent equations which 
derive from (|37|) for this geometry. This is valid for arbitrary A(r) = 1 + 
a{r) and B(r) = 1 + b(r). In the MOND regime one has \a(r)\ 1 and 
\b(r)\ 1, for which case we show that the equations can be simplified 
without making any assumption about the interpolating function T(x). In 
terms of the general weak field ansatz (J20|) . these simplified equations prove 
that ei + €2 = 0. It follows from section 2 that no model of the form ((7|) 
can be consistent with galaxy lensing. However, it is possible to choose the 
interpolating function jF(x) so as to reproduce MOND's success with galactic 
rotation curves. The section closes by doing this. 

The general spherically symmetric and static geometry (jl2H4|) gives rise 
the following Ricci scalar, 

B" B' (A' B'\ 2 (A' B'\ 2 / 1\ 

r = -a^ + 2Ab{-a + 13) + 7a{-a--b) + vA 1 'a) ■ (44) 
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In this geometry, and acting upon a function only of r, the covariant d'Alem- 
bertian reduces to, 



d f 2 IB d 



" 2 \fAB~dr\ V Adr 



(45) 



The differential equation which defines the small potential therefore takes 
the form, 



> \B , 



aKa + b)- 

(46) 

Assuming the two parenthesized terms vanish at r = we can write, 



^^^(^-l^VV^-i)^!) . (47) 
The differential equation that defines the large potential is, 



Assuming again that the parenthesized terms vanish at r = we can write 



(48) 
ite, 

(49) 



And those of its second covariant derivatives we shall need are, 



B' 

$ tt = 7$' , $T 



4' 1 



.4 



_ 2_ IfB' A' 
r 2\B A 



■ (50) 



In this geometry only the diagonal components of the field equations (|37j) 
are nontrivial. The 99 and 00 equations are proportional to one another, 
and conservation gives both from the tt and rr equations the same as for the 
stress-energy, 



Ta, 



sin 2 (#) 



r 3 B l A 

Tee ~ 2A\2BB Ttt + 



d 2 A' B' 
dr r A 2B\ 



T 



(51) 
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The system can therefore be defined by its tt (times A/B) and rr equations, 



RirH A A A n 2 A' 

^T tt = 2$" + U' + gG tt (l - 2$) + ^AT - ±# - J* , (52) 



87TG, 



4 



B' 



+ G rr (l - 2$) - AT - 
r 2c 4 £> 



The tt and rr components of the Einstein tensor are, 



A G - A ' + 
S Gtt " rA + 



A 



B' 

G rr — 
rB 



A - 1 



(53) 



(54) 



We now evaluate the derivative of the small potential (|47|) to leading order 
in the weak fields, a(r) and b(r), 



r 



(55) 



One might worry that the integral in (|47|) contributes as well, but note that 
the integrand exactly vanishes for A = B~ x . This means that the integral 
cannot contribute much in the regime for which general relativity applies. 
The integrand is nonzero in the MOND regime, but it is also second order 
in the weak fields, a and h. We can therefore ignore this term altogether. 

In the asymptotic regime we can assume that each derivative adds a factor 
of 1/r. Hence <p'(r) goes like 1/r times the small numbers a(r) or b(r). It 
follows that $'/r is much larger in magnitude than <//$'. By similar reasoning 
we recognize that $'/r and $" dominate the other MOND corrections, 



> 



~B 



(56) 



In the weak field limit it seems reasonable to assume T rr = while still 
allowing a nonzero A/BT U = p. We still don't know how the two leading 
MOND terms (|56J) on the gravitational side of the equations compare with 
the leading terms from general relativity. Including both gives, 



4^. a' a 



a 



+ ... 



8ttG 
. 



p(r) 



(57) 
(58) 
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The first of these equations (joT|) can be integrated to give, 

+ 72 + • • • = + ^I-T / *"r' 2 p(r') • (59) 
Adding (j3Hj) and (JoTJj) cancels the leading MOND corrections, 

- + - + ■•■ = - + "XT" / d^pif) . 60) 

^? gal 

Equation (jHOj) is interesting because it has no dependence upon the still 
unknown interpolating function Fix). We can therefore use it to make state- 
ments about all models of the type (J7J). Under the assumption of no dark 
matter halos, the mass integral must eventually stop growing, in which case 
the left hand side falls off like 1/r 3 . So if b'(r) goes like a constant times 
1/r then a(r) must go like minus the same constant. In terms of the generic 
ansatz (|20|) of section 2 we have just demonstrated that all models of the 
type (J7J) have e\ + £2 = 0. As discussed in section 2, this means that galaxies 
without halos of dark matter would give far too little lensing. These models 
do not give phenomeno logically viable realizations of MOND. 

It is still interesting to see if the interpolating function T[x) can be chosen 
to reproduce MOND rotation curves. For this purpose let us consider a sphere 
of radius R with very low, constant density, 

3 Mr 2 

p{r) = ^ 9{R - r) - (61) 

If the density is small enough the MOND regime prevails throughout, as in 
a low surface brightness galaxy. This means that (JH7|) can be integrated all 
the way down to r = to give, 

2$' + - + ■■• = ^ fdr'r' 2 P {r') . (62) 
r c 4 r z J 

We can also use (J55|) to eliminate b'(r) in — r times 

W + (f/ -- + ... = . (63) 
r 

Now eliminate a(r) by adding (jffiZj) and ([63)1 . and then use (|49jl to obtain an 
equation for the small potential, 

c 4 r 5 



«0 



+ • • • = / dr'r' 2 P (r') . (64) 
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For r > R the mass integral is constant, 

- 4 <-' 2 - 2GM 



Now recall from section 2 that MOND requires e 2 = 2, and we have just seen 
that any model of the class (J7J) must have ex = — e 2 . The weak field limit 
(155(1 for the small potential therefore implies we must have, 



6 a GM 

( f( r )^--\l^^ + --- ( 66 ) 

It follows that the constant term within the square brackets of must 
exactly cancel, and that the next order term must involve one power of ip'. 
Working out the algebra gives, 



T ' {x) = - 5 " is + ° w =* = "i " m + 0{x2) ■ (67) 

The associated weak fields are, 



AGM a GM 8GM a GM f r\ 



For the general weak field ansatz (j2"Uj) of section 2 we have just shown —25 



8 2 = -§ and — ei = e 2 = 2. 



A potentially troublesome point is that equation (J55j) involves terms of 
second order in the weak fields a(r) and 6(r), 



, ^/cVm c 2 . 2 2GM , , 

\ a /J 18an err 2 



But in (jHEJ) we previously neglected such second order terms to derive the 
simplified equations (J57H58j) which pertain in the MOND regime. Closer 
inspection reveals that all of the terms neglected in (|56|) contribute terms to 
the left hand side of ()69|) which are small for r <C i?hor ~ 10 26 m, 

c 2 

rw' 2 < — y/ 2 ~ 10 27 m x in' 2 . (70) 
a 

Note also that one must involve quadratic terms like those of (|59*jl in order to 
make the weak fields go like the square root of the system's mass, as MOND 
predicts. 
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Enforcing the MOND limit determines only the first two terms in the 
small x expansion of the interpolating function T{x). There are many ways 
of extending this to a formula for general x. The only phenomenological 
constraint is that we need the MOND corrections to be acceptably small 
in the general relativistic regime of large x. For example, we can make 
T(x) — ► — y|x|2 for large \x\ with the following extension, 

= -igw + Jgw , (71) 

1 + eW 5 U + \\x\i\ 

22 I I — 7 I I 

F{x) = _ TFl 2 + g\ x \ + Uln(l + -\xfi) . (72) 
l + ||x|2 V 6 / 

For |x| 3> 1 this would typically suppress MOND corrections by some char- 
acteristic length of the system divided by c 2 /ao ~ 10 27 m. If that is not 
sufficient one can always extend J-{x) differently to obtain more suppression. 



5 Homogeneous and isotropic sources 

Although our relativistic formulation of MOND is ruled out by lensing it 
seems a pity not to work out the cosmology now that we have the formalism. 
The exercise also affords a potentially important caveat on just how much a 
general formulation of MOND can change in passing from the static geome- 
tries of galaxies to the time dependent geometry of cosmology. We begin by 
working out the small and large potentials for a homogeneous, isotropic and 
spatially flat metric, 

ds 2 = -c 2 dt 2 + a 2 {t)dx ■ dx . (73) 

The nonlocal field equations ()37|) are next specialized to this geometry. Then 
the MOND limit is taken. The section closes by considering what happens 
when matter domination follows a long period of radiation domination. 
In this geometry the Ricci scalar is, 

H = -. (74) 



R(t) . (75) 



c 2 R = 6H+ 12H 2 where 
The small potential is defined by the equation, 

Dip® = -a-^(a^) = 
Yy ' dct\ dctJ 
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If we define the initial values of ip and its first derivative to be zero, it takes 
the simple form, 

(p(t) = - J" dt'a-'\t') J* dt"a\t") (6H{t") + 12# 2 (i")) . (76) 
The large potential is defined by the differential equation, 

d*(j=9Sr*r) = d»(V=99^F') =► j^) = J t {^') • (77) 
If we again assume null initial values the result is, 

= j^dt'^t')r{-c 2 af^ 2 {t')^ . (78) 

The nonzero components of the second covariant derivative are, 

$ ;00 = cr 2 $ , = -c- 2 H$ 9ij . (79) 

We assume the stress-energy tensor to take the perfect fluid form, 

T nv = P9^u + (p + p)u^u v . (80) 

Stress-energy conservation implies, p = —3H(p + p). The nonzero compo- 
nents of the Einstein tensor are, 

c 2 G 00 = 3H 2 , (81) 
c 2 G il = ~(2H + 3H 2 )g ij . (82) 

We therefore extract two nontrivial relations from the general equation (|37|). 

n 2 

8nGc- 2 p = -6#$ + 3H 2 (1 - 2$) + , (83) 

2c 2 

•■ -a 2 

8nGc- 2 p = 2$ + AH<& -(2H + 3H 2 )(l- 2$) - 0$ \T ■ (84) 

As in general relativity, only one of these equations is independent. The 
second equation ()84|) follows from the first and stress-energy conservation. 

In the MOND regime the interpolating function and its derivative take 
the forms, 

1 : 1 

F{x) — ► , T{x) — ► --sgn(x) . (85) 

6 6 
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For cosmology the argument x = — (cip/ao) 2 is negative so the large potential 
has the same sign as the small potential, 



Mt) ■ 



(86) 



In the MOND regime we can therefore express (jHEj) as, 

-H^ + 3H 2 (l- -tp) - ^p 2 + ■ ■ • = 8nGc 



P 



(87) 



An additional specialization of great interest to cosmology is the case of 
a power law scale factor, 



a{t) = (l+Hity 



(88) 



Here Hi is 1/s times the Hubble parameter at t — 0. Substituting into (f?7)|) 
gives the small potential, 



<p(t) = —6s 



2s 



3s -1 



In 



1 + Ht 



'I - 3s 



1 + Ht 



l-3s 



1 



(89) 



At late times only the logarithm matters. In this regime we can also express 
p in terms of the Hubble parameter, 

'2s - 



<p(t) 



-6 



3s 



H(t). 



(90) 



We can therefore write the MOND analog of the Friedman equation for power 
law expansion, 



3 l + 2ff - a 2 + 2saln 



1 + Ht }H 2 (t) + ■■■ = 8nGc- 2 p(t) 



(91) 



where a = (2s — l)/(3s — 1). The effect of the logarithm on s > | power 
laws is to gradually slow the expansion. This makes rough physical sense if 
we think in terms of MOND strengthening the force of gravity in the weak 
field regime. 

For the case of radiation domination (s = 1/2 and a = 0) we note that 
(p(t) = 0! The large potential also vanishes — exactly, not just in the MOND 
regime. Hence the equations reduce to those of general relativity, but with 
the energy and pressure coming from only ordinary matter. This is phe- 
nomenologically unacceptable. One of many things that goes wrong is nu- 
cleosynthesis. 
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6 Discussion 



We have succeeded in embedding the MOND force law in a set of covariant, 
causal and conserved field equations for the metric. The model suffers from 
at least two fatal phenomenological problems in which its predictions agree 
with those of general relativity without dark matter. The deflection of light 
and cosmology during radiation domination both have this property. Of 
course these problems do not necessarily mean that MOND is wrong, only 
that our realization of it is. 

Our model consists of corrections which are based on the small potential 
ip[g] = D -1 /?. One naturally wonders if it is possible to find a nonlocal scalar 
potential that avoids the problems with lensing and cosmology while still 
keeping the MOND force law. For example, one might replace the covariant 
d'Alembertian with the conformal one, 

<Pc[g] = ^r-R where n e = U-^R. (92) 

The distinction between □ and C c disappears in the weak field regime be- 
cause R goes like one power of the weak fields over r 2 . So such a model would 
still give the MOND force law. Unfortunately it would also give too little 
lensing precisely because its weak field limit agrees with that of <p[g]. Be- 
cause ip c vanishes with R, it would also have problems dealing with radiation 
domination without dark matter. 

The fact that R vanishes for a radiation dominated universe means that 
we should avoid it as the source upon which the nonlocal operator acts. The 
next most complicated scalar potential would seem to be, 

<P2[g] = ~(R>"R»J) . (93) 

Because (p% has roughly two derivatives acting upon two powers of the weak 
fields, one must also change the Lagrangian, 



£ 2 



c 4 



16nG 



R + c- 4 al^2[g})]V^g ■ (94) 



The interpolating function JF 2 (x) would become linear in the MOND regime. 

One interesting thing about our model is that it becomes conformally 
invariant in the MOND limit. The first hint of this came when all our MOND 
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corrections cancelled out of the formula for the deflection of light. Photons 
are conformally invariant so they experience no deflection due to conformal 
transformations of the metric. 

To prove asymptotic conformal invariance, note that in the MOND limit 
{F{x) — ► — qX for x > 0) the large potential is, 

*\9\ = ^{v*F')^-\<P- (95) 
In this limit the field equations ()37|) take the form, 

S-kGc^T^ = - (g^U(p - (p.^ + Rp, v - -^g^vR + ■■■ ■ (96) 

The right hand side is traceless, which means that the linearized theory is 
conformally invariant. 

This sort of asymptotic conformal invariance was also found by Bekenstein 
and Milgrom ^SJ|2Z1- Conformal invariance was certainly not built into either 
model. For example, the trace of our full field equations (J37j) is nonzero, 

2a 2 

8nGc- 4 T^ = -6n$ - R[l - 2$] + 2<p'»<$> tfl + f^f^J 7 ' - -^T . (97) 

Asymptotic conformal invariance arises from enforcing the MOND limit. 
Since MOND requires the weak fields to go as the square root of the source 
mass, it is necessary that terms linear in the weak fields drop out of some 
component of the field equations. The only distinguished component in a 
conserved, tensor equation is the trace. Hence terms linear in the weak fields 
must drop out of the trace of the field equations, which means that the 
linearized theory is conformally invariant. 

If asymptotic conformal invariance is generic it means that no metric- 
based formulation of MOND can give the required amount of gravitational 
lensing. It might also bear on the view that the successes of MOND derive 
from galaxy formation and evolution flowing, through conventional physics, 
towards some sort of fixed point. This is because critical phenomena and 
universality are typically characterized by conformal invariance. It might be 
interesting if the analogy could be pursued sufficiently to generate quantita- 
tive predictions. 
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